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Pee OUUCT TON 


m@eemotion of bodies through fluids and fluid motion in 
Semebout bodies gave rise to a great number of practical and 
theoretical studies regarding the determination of resistance 
in steady and time-dependent flows. 

Jbservations as well as numerical experiments show that 
the wake of a bluff body is comprised of an alternating vortex 
street. The character of the vortices immediately behind the 
body and in the wake further downstream depend, for a steady 
ambient flow, on the Reynolds number and the intensity and 
scale of the turbulence present in the ambient flow. For a 
time-dependent flow, the instantaneous state as well as the 
past history of the flow play significant roles and it is not 
possible to give a general set of normalized parameters. 
meets tic reasoning, and laboratory and numerical experiments 
Seem Will have to supplement the information gathered from 
melatively idealized solutions in order to correctly identify 
the most important governing parameters. 

mvectets the £lWid is assumed to be inviscid, flow involv- 
ing separation is beyond the reach of rigorous calculation. 
Im reality, the flow behind a bluff body moving steadily 
menomen a tluid (or held stationary in a steadily moving 
fluid) is accompanied by large scale unsteadiness. Further- 


aeomeany type Of disturbance or unsteadiness experienced 


JES) 





by the ambient flow and/or the body introduces additional 
mmees In the characteristics of the flow and in the fluid- 
mermecure interaction. 

The formation of a wake gives rise not only to a form 
drag, as would be the case if the motion were steady, but 
meso tO Slgnificant changes in the inertial forces. The 
velocity-dependent form drag is not the same as that for 
mmeedy flow in a viscous fluid; and, the acceleration- 
@ependent inertial resistance is not the same as for un- 
eeearaced flow in an inviscid fluid. In other words, the 
drag and inertial forces are interdependent as well as time- 
dependent. Thus the relationships between the various re- 
Sistance components must be determined in terms of the 
unsteadiness of the ambient flow. Also, these relationships 
femeeake into consideration the geometry of the body, the 
@eearee Of upstream turbulence, the roughness of the body, 
Seaethe past history of the flow. 

The hydrodynamics of unsteady flow, in particular those 
Set 1n motion impulsively from rest, has been most aptly 
described by Sedov [1]. Birkhoff and Zarantonello [2] have 
Drought together results of nearly a hundred years of research 
On the motion of jets, wakes and cavities. More recently, 
Mevouutionary changes in computers and numerical methods have 
generated new concepts which are scattered throughout the 


literature. 
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As in a great many steady flow problems, where potential 
flow analysis provides many approximate solutions, real wake 
behavior using models involving potential theory have been 
tried. (For example, see Sarpkaya and Shoaff [3].) Some of 
these approximations may be expected to provide realistic 
Solutions and make it possible to eventually introduce the 
effects of roughness, viscosity, gravity, turbulence and past 
history, among other variables. 

Favorable comparison between theoretical results and ex- 
perimental observations may not always occur, but the objec- 
tive is not always the achievement of a completely favorable 
@emiparison. Rather, it is the careful selection of those 
theoretical predictions found to be reasonably accurate (in 
their agreement with experiments) that should be regarded as 
meemrst Step toward further exploration of potential flow 
methods and then eventually the complete modeling of all 
phases of the flow phenomena. These efforts logically lead 
mothe reasons why impulsive flows are studied extensively. 

Bemmoted earlier, time-dependent flows involve changes 
MieooOth the drag and inertial forces. The history of the 
Merton may affect significantly the instantaneous value of 
the resistance and the time it takes to reach a given value. 
Thus, a general time-dependent flow cannot yet be analyzed 
mamerically and, therefore, does not provide an opportunity 
for the development and testing of numerical techniques and 


for the comparison of the results with those obtained 
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meperimentally. However, the impulsively-started flow about 
Secample bluff body is relatively more manageable and can be 
Smalyzed analytically and/or numerically, at least for the 
Gummyeotages Of motion, and provides an excellent opportunity 
eemene Validation of the methods used through the comparison 
Seememtneoretical and experimental results. [t 1s partly 
meeeise Of this reason that the study of impulsively-started 
flows has attracted a great deal of attention during the past 
Mmety years. 

meoreestarted impulsively from rest are of importance from 
meee tical point of view also. Examples of applications of 
Meee }OWS include the effects upon bodies immersed in a 
meme water reactor during a loss-of-cooling-water accident 
Seamene understanding of the motion of missiles at high angles 
Saeattack. As will be noted later, impulsive flow is analagous 
Bemene evolution of separated flow about slender bodies moving 
Meeaeh angles of attack. Underwater vehicles, such as torpe- 
does and submarines, have impulsive changes imposed upon their 
Meron, particularly during a turn or dive. 

Mastly, from a basic research point of view, the evolution 
Saecme wake and the formation of the first few vortices are 
far more important than the subsequent steady state. As the 
mesults presented herein show, the most dramatic changes in 
memrOorce-transfer coefficients take place during the early 


Stages of motion. Understanding of the kinematics and dynamics 
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Bemtene carly stages of the evolution of the wake provides 
mean le information regarding the growth of shear layers 
meemene distribution of vorticity. 

fers an view of its practical and theoretical signif- 
meee that this investigation of impulsively-started flow 
Meemeechiree types of bluff bodies (a D-shaped cylinder, a 


Mmeemeped cylinder and a flat plate) has been undertaken. 
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eee Ore lie PREVIOUS INVESTIGATIONS 


Impulsively-started flow 1s one of the unsteady flow situ- 
ations where analytical and numerical solutions exist at least 
for small times and relatively low Reynolds numbers. 

Mune the early stages of motion the vorticity does not 
have enough time to diffuse. The boundary layers are very 
thin, the flow is essentially irrotational and the fluid force 
meme On the body 1S primarily inertial. The inertia coef- 
Meetent is Cm = 1 + Ca where Ca is the added mass coefficient 
Seeained from potential flow theory. For bodies without sharp 
Mutem>s, aS in the case of a circular cylinder, separation 
Goes not occur immediately. Furthermore, the separation does 
not necessarily initiate at the downstream stagnation point 
(as in the case of an elliptic cylinder). 

For two-dimensional cylinders, it can be shown that the 
mem@eration begins, after a time t* at a place where the abso- 
lute value of dU/dx is largest. The relationship between t* 
ema Du/dt is [4] 


1+ (1 + wo) Pit = 0 (1) 


meee smooth circular cylinder started impulsively from rest to 

Meeonstant velocity, the distance covered until separation begins 
Mees = 0.3551C (C is the radius of the cylinder). Here, the sep- 
aration does begin at the rear stagnation point. For a uniformly 


meeeteratins circular cylinder the corresponding distance 


Meaveled is s = 0.52C. 
20 





FOr axisymmetric bodies t* is given by the expression [4] 


dU 4 Wied. 
* 2==—<—= = —_ 
ees | dx oe aT 7 Oh Noes dx 0 (2) 
Seemaespnere impulsively set in motion s = 0.392C. The dis- 


tance covered by the sphere until the onset of separation is 
somew.at larger, as inthe case of the cylinder, when the 
Meeere iS accelerated uniformly from rest. The rate of 
acceleration as well as the history of acceleration is 
mgorecant in the calculation of the relative distance 
Sovered prior to the occurrence of separation. For bodies 
With sharp corners separation starts immediately. 

Consider an impulsive change superimposed on an already 
established flow pattern and, with that change, the role played 
by separation on the added mass. Just prior to the impulsive 
Change the drag coefficient is given by its steady state value 
at the corresponding Reynolds number. Sears, as reported by 
Rott [5], has shown that "the initial motion following the 
impulsive change of the conditions consists of the super- 
position of the velocity pattern existing just before the 
Change and the inviscid flow velocity pattern due to the 
impulsive boundary values.'' In other words, at the initial 
instants of impulsive change the drag coefficient is equal 
Beetts steady state value and Cm = 1 + Ca (Ca again comes 
meom the potential theory; for example, Cm = 2 for a Cir- 


Baear cylinder). As time progresses, neither Cd nor Ca 
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remains the same and both change with the evolution of the 
meow ihese changes are affected by the past history of the 
motion and the features of the current state. Thus, the 
changes in the added mass coefficient come about not only 
because of separation but also because of changes in the 
state of the separated flow, such as the additional move- 
ment of separation points and the increase or decrease of 
M@emracte Of circulation. 

Theoretical investigations of impulsSively-started motion 
Seeeeecircular cylinder in a fluid otherwise at rest were 
feet conducted in the early 1900's and confined to very 
small Reynolds numbers. Such a motion was first considered 
by Blasius [6] in 1908 and his work was later extended by 
Goldstein and Rosenhead [7] in 1936, Gortler [8,9] in 1944 
merges, Schuh [10] in 1953, Watson [11] in 1955, and Wundt 
Meeiein 1955. It was found, as noted earlier, that after a 
@emtain lapse of time, the boundary layer separates from the 
Surface of the cylinder; the time and location of separation 
depending on the Reynolds number and the bluffness of the 
body. The separation points then move rapidly around the 
cylinder until at some greater time they coincide with the 
average positions of the points of laminar Separation for 
Steady flow. 

Finite difference techniques have been employed in the 
analysis of separated flows by several investigators (Payne 


[13] in 1958, Hirota and Miyakoda [14] in 1965, Kawaguti 
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meemeain [15] in 1966, Wang [16] in 1967, Jain and Rao [17] 
fmeoo9, Rimon [18] in 1969, Son and Hanratty [19] in 1969, 
Seeman and Szewczyk [20] in 1969, Honji [21] in 1972, Mehta 
meemeeenan |22] in 1972, Collins and Dennis [23,24] in 1973, 
Seemed thompson [25] in 1973, Lugt and Haussling [26] in 1974, 
Memonis and Tsahalis [27] in 1974, Bar-lev and Yang [28] in 
ioe Panniker and Lavan [29] in 1975, and Cebeci [30] in 
ime). Among others, Tuann and Olson [31] employed the finite 
element method in 1976. 

Experiments at relatively low Reynolds numbers have been 
Memerced by Schwabe [32] in 1935, Taneda and Honji [33] in 
1972, Taneda [34] in 1972, and Coutanceau and Bouard [35] in 
77. 

Very little experimental data for impulsively-started flow 
at sufficiently high Reynolds numbers (for Reynolds numbers in 
Mpmoupercritical and post-critical regimes) exists. This is 
partly because of the experimental difficulties encountered in 
establishing a vibration-free impulsively-started steady flow 
Emempeartly because of the instrumentation required to measure 
Maenwtransient quantities involved. In fact, the force that 
acts on a cylinder in impulsive flow at relatively high 
Reynolds numbers has been measured directly only by Sarpkaya 
mon 57| in 1968 and 1978. 

Schwabe [32] used a circular cylinder with a radius of 
© - 1.7 inches, and conducted his experiments in an open 


Water channel. The velocity of the cylinder was U = 0.528 
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anches per second at a Reynolds number of about 600. Veloci- 
Meesewere determined from the path lengths of particles sus- 
Menmded On the water surface. Other quantities such as the 
time rate of change of the difference of the square of the 
meemeres, rate of change of the circulation, and the radius 
Mumpetliyvacire Of the streamlinesneeded to calculate the drag 
coefricient from the Bernoulli equation were all determined 
Meeeerapiical methods. A careful examination of this work 
mu@meactes that a considerable amount of experimental error 
May have existed in the evaluation of the pressures and, 
@meeetore, in the resulting drag coefficient. (Schwabe's 
Grae coefficient is about twice the steady-state value and 
meme increaSing, when the cylinder has moved about 9 body 
madii.) 

Bangham et al. [38] carried out a number of experiments 
Baea Shock mb to observe the influence of Reynolds and Mach 
memoeer On the impulsive loading of a 0.5 inch cylinder for 
Reynolds numbers in the range of 31,000 and 77,000 and the 
Seeeeiumoecr range of 0.15 to 0.4. The pressures were deter- 
mined from the density fields and the drag coefficients were 
Mmeetlated by an integration of the pressures. The variation 
Meeettie drag coefficient with UT/C, for a given set of Re and 
Ma values and from one set of Re, Ma to another, is signif- 
Meant. These variations occur partly because of the growth 


meaesnedding of vortices and partly because of different 





rates of decay of the initially high pressure on the rear of 
the Cylinder due to passage of the two branches of the orig- 
mia shock. 

Friberg [59] quickly immersed a circular cylinder into 
the steady and uniform free-stream flow field of a water 
table. Asher and Dosanjh [40] conducted experiments in a 
shock tube in a manner similar to that done by Bingham et al., 
and c2termined the characteristics of the wake. (For example, 
the position of the vortices, their relative velocities and 
Strouhal number). Their experiments suffered from the same 
drawbacks as those of Bingham et al. Other than through 
mumeérical simulation, there is no mechanical or traveling 
mmeek system which is capable of generating a truly impul- 
Mme low. Efforts to create impulsive or uniformly accel- 
erated flow at high Reynolds numbers or accelerations in a 
Brquid medium may be hampered by the generation of compres- 
Sion and rarefaction waves and regions of intense cavitation. 
These have been some of the difficulties of experiments with 
impu_sive flows. 

Sarpkaya [37] recently repeated his previous work with 
Serecilar cylinders using a water tunnel with considerably 
more-sophisticated instrumentation, and simultaneously 
measured both the drag and lift forces. Sarpkaya and Shoaff 
[3,4] attempted to numerically simulate the evolution of the 
memery Stages of motion through the use of the discrete vortex 


model. 
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Mmeey have shown that the drag coefficient in the initial 
merees (UI/C = 4) of an impulsively-started flow can exceed 
imeemsceady value by as much as 30 percent. In the early 
memmeas Of the flow, vorticity is slow to diffuse and, there- 
Mememeaccumulates rapidly in the close vicinity of the cyl- 
inder. Although the growing vortex soon reaches unstable 
Meeperrions and separates from its shear layer, the growth 
Or tne vortices is so rapid that they become larger than 
their quasi-steady-state size (at about the same positions} 
Derore they separate from their shear layers. This leads to 
miemtarce drag coefficient observed in the experiments (see 
Roos and Willmarth [42] for similar observations with spheres). 
Summeeey atter the onset of asymmetry, the drag coefficient 
Memeeeases Sharply and the lift coefficient begins to increase. 
Seeeequently, the lift coefficient oscillates with the fre- 
Gmency of the MecdcimommoOneLeces, from the same side of the 
eylinder. 

maeime Subsonic to moderately supersonic velocity range, 
miemapproximate flow similarity between the development of 
Mm@encross~flow with distance along an inclined body of uni- 
meme diameter and the development with time of impulsive 
flow about a cylinder is known as the "cross-flow analogy." 
This analogy was first suggested by Allen and Perkins [45] 
in 1951. It has been subsequently used by many other 
Meeecarchers to calculate the in-plane normal force and the 


mme-or-plane force (side force normal to the plane of flight) 
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mere On Slender bodies moving at high angles of attack. A 
detailed discussion of the analogy, along with extensive meas- 
urements for various nose shapes and body combinations, may 
be found in Thomson and Morrison [44,45]; Thomson [46]; Bostock 
Seeeeeveamont [48]; Lamont and Hunt [49]; Wardlaw [50]; and 
memesson and Reding [51]. 

\ blunt-nosed cylinder, such as a submarine, at high angles 
Of attack generates a stationary asymmetric vortex array which 
meeeieiar to the Karman vortex street. Here, the approximate 
mies pace equivalence for the analogy is possible because the 
vortices have an axial degree of freedom for their lift-off. 
memever, it must be pointed out that the analogy is far from 
Memerect. On a pointed slender body the first asymmetric 
Meeeex Palr Originates at the apex and the analogy no longer 
mpplies. Asymmetric vortices give rise to a coupling between 
longitudinal and lateral degrees of freedom which is often 
@iscontinuous and associated with hysteresis effects. Con- 
Mementliy, the phenomenon ceases to be a simple two-dimensional 
Space-time equivalence and becomes a complex three-dimensional 
muedestructure interaction problem. Additional discussion 
@umeceierences for this and other types of time-dependent 
flows are presented by Sarpkaya and Isaacson [52]. 

Because most of the past analytical and experimental 
@eudtes dealt with a circular cylinder attention herein has 
been, thus far, primarily directed to impulsive flow about 


that particular body. Very few studies exist regarding other 
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@meeoodies. Fage and Johansen [53], in 1927 and 1928, con- 
@ucted experiments with various bluff bodies, including the 
meeeaeion Of flow about an inclined flat plate. From their 
mieenious efforts it is known that vorticity is shed from the 
two sides of an asymmetric body at the same rate and, that 
mieemotioOn in a vortex sheet is steady near the body, except 
possibly near the inner edge of the shear layers. They have 
also provided detailed observations regarding the growth and 
mepemsion of the shear layers. 

In 1975, Sarpkaya [54], presented an inviscid flow model 
Of two-dimensional vortex shedding for an impulsively-started 
flow over an inclined plate. This work was followed by Fink 
mm@mooh {55| with a slightly refined version of the discrete 
vortex model. Their analysis was confined to the early stages 
Mmemene vortex development and the calculations were not ex- 
mended for sufficiently long times to examine the shedding 
of the vortices from the plate. Belotserkovski and Nisht [56] 
presented the results of a discrete vortex model of flow over 
a flat plate at large angles of attack. They did not describe 
Mmemcetails of their analysis. Their calculations resulted 
memomrealistically large oscillations of the drag coefficient. 
Mmewcarly stages of the evolution of vortices from a flat 
plate was studied by Wedemeyer [57]. Subsequently, Pullin 
[58] extended this work to flow past wedges at various in- 
Cluded angles. Telste and Lugt [59] examined the vortex 


Shedding from finned circular cylinders and found that the 
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memee coeificients are very sensitive to the numerical model. 
memever, they have not reported those force coefficients. Kiya 
and Arie [60] attempted to extend Sarpkaya's work [54] by assum- 
Mmochat the point vortices may be introduced at fixed points 
Meee tips of the plate. Their calculations resulted in 
feeeeestiuctuations in the drag coefficient. This and the con- 
Meeeure that the nascent vortices are shed from fixed positions 
Mmevents realistic assessment of their results. 

memmoted before, numerous investigations have used the 
Navier-Stokes equations, together with appropriate finite- 
difference methods, and modeled the evolution of the wake 
Guring the earliest stages of the motion. Unfortunately, 
such calculations are limited to very small Reynolds numbers 
(in the order of 100); partly because of problems encountered 
with numerical Meatless wand partly be@ausesof the size 
Meittations of the existing computers. Consequently, the 
Mesmuits of these studies cannot be compared with those ob- 
tained experimentally at larger and practically more signif- 
icant Reynolds numbers. 

the work described herein has been conducted for the pur- 
Meese Of providing data at significant Reynolds numbers for a 
Meme plate and for two cylindrical shapes which have not been 
analyzed or tested previously. [It is hoped that additional 
numerical work will be undertaken in the future and will be 


guided by the results presented herein. 





Pie ec rER LueNn LAL SD EGU PMENT wiNDe PROCEDURES 


Bee VERTICAL WATER TUNNEL 

m@emexperiments were conducted in a vertical water tunnel 
(see Figure 1) which was previously used for the study of im- 
mebemvertlOW about circular cylinders in 1977 [61]. A quick- 
release valve located at the base of the tunnel is used to 
@peate an impulsively-started flow of constant linear velocity. 
A partial drawing of the mushroom-like seating surface of the 
quick release valve is shown in Figure 2. 

mmoracr tO prevent distortion of the force and moment in- 
M@meations, the water-side profile of the mushroom valve has 
been especially designed to ensure continuous undisturbed flow 
Meeeencne Seat while the valve is in any open position. As 
Meme im Figure 2, when the valve is closed, it is in the 
mely Up position. It seats against an ‘'0' ring inserted on 
MemDOttTOm Of the seating surface so that no leakage is pres- 
eeorior tO initiating fluid motion. 

The vertical position of the mushroom valve is controlled 
by a three-way valve mounted beneath the tunnel (Figure 5). 
The stem extends downward from the mushroom valve and is 
directly coupled to the control valve piston assembly. Com- 
Meeosed air is provided to the two air chambers in the upper 
part of the valve. A two-way valve in the air supply line 


meovides on-off control of the quick release valve and, 
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Two-Way Rir Supply Valve. 


. GQuick-Release Valve (inside the Base). 


2 
3. Moment Indicating System. 
4, Urag7Lift Force Indtcating System. 


5S. Front Recess and Shadow sox. 


Figure i. Vertical Water Tunnel. 
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@heretore, the flow itself. A Sight glass provides visual 
Memeeoring of the fluid level during an experiment. 

Ween Opening the two-way air supply valve, the differen- 
mereepressure between the two air chambers in the upper part 
Meee control valve initiates motion of the piston, rapidly 
Meeming the mushroom valve. Thus, flow is initiated and the 
Sumancorains Out into a reservoir beneath the tunnel. Subse- 
@m@ent valve motion is regulated by the vertical motion of the 
Meeeeomein the lower part of the control valve and the viscos- 
MueeoreOll in the liquid chamber. The area of the opening 
(and consequently the amount of resistance which the piston 
Smeounters) between the liquid chamber and the upper air 
Siamber can be varied by opening or closing the dual ports 
moeene piston (Figure 4). Oil viscosity and supply air pres- 
Sure can also be adjusted. 

These adjustments allow constant velocities at desired 
Mmeees tO be obtained. Following the rapid initial opening, 
which accelerates the flow in 0.1 seconds or less, slower 
Surtner Opening of the mushroom valve sustains a controlled 
drop of the tunnel water level. This controlled slowdown of 
Meeerate at which the mushroom valve opens, provided steady 
memeeevelocities of about 0.9 feet per second for the experi- 
Ments described herein. No attempt was made to repeat the 


Meperiments at other velocities. 
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B. VELOCITY, FORCE AND MOMENT MEASUREMENTS 

Velocity was measured through the use of a variable resist- 
ance probe. A five foot long platinum wire, placed vertically 
in the tunnel and mounted away from the walls, provides tunnel 
Water level indication to an amplifier-recorder assembly. 
meror tO Conducting any experiments, impulsive flow was ini- 
tiated several times. Adjustments to the quick release valve 
control system were made, as necessary, to ensure a linear 
meeme On the elevation versus time plot on the recorder. 
Welocity calibration settings were thereby verified. Addi- 
Mumiy, fOr each run, a Consistent water level of three 
Beet above the lower end of the platinum wire was visually 
meemising the installed sight glass. This enabled run to 
@emmeconsistency of the velocity profile to be verified. 
meeure 5 shows an example of the recorder output for velocity 
Meerenent. 

Peeeelectronics 10 ke load transducers were used to mea- 
mires tne instantaneous lift and drag forces and moment on the 
Bedies tested. Special housings were built for each gage so 
meat they could be mounted on the tunnel at each end of the 
best body. As shown in the tunnel photograph (Figure 1), 
the moment was measured from one end of the test body while 
the lift or drag was simultaneously measured at the other end. 
The bellows protecting the strain gages were filled with Dow- 
Beriing RIV coating for water proofing and then the ends of 


the bellows were sealed air tight. 
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Figure 5. Analog Record for the D-Shaped Body at d deg. 
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Memordcer to measure drag or lift without affecting the 
moment measurement, the test bodies were mounted in a self- 
aligning bearing (Figure 6) whose housing was connected to 
Mm@eetorce transducer. This allowed the test body to be 
mreely mounted (thereby permitting true moment indication 
@t the other end) and accurately aligned. The bearing and 
housing for the drag and lift gage assembly were also de- 
Signed to let the gage be rotated to either measure drag or 
mrt . 

Peeter mOunting the test body, the force gage lift posi- 
ti0n Was determined by rotating the strain gage beam to the 
Wertical position and then hanging different loads from the 
memgye ihe lift position was fixed when a zero drag force 
reading was obtained. The final position was marked with an 
metenment hole, drilled through the indicating plate and 
housing at the outer end of the assembly, at the top of the 
Wertical axis of the indicating plate (Figure 7). A set pin 
@emeeeten Placed in this hole to lock the gage in the lift 
Mesttion. Another alignment hole, marking the drag position, 
Was drilled into the indicating plate exactly 90 degrees from 
Meeevertical. Four bolts in the gage housing hold it rigidly 
in place. Removal of these bolts and the set pin permitted 
POtation of the gage to either the drag or lift position. 
Mien the pin is placed in the proper location to secure the 


alignment accurately. 
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pe the moment side of the tunnel, the test body is con- 
nected to another 10 kg load transducer. However, in this 
case the load transducer acts as a connection to the strain 
gage (torque-beam) which actually measures and transmits the 
moment indications. As displayed in Figure 8, the water side 
(front) of the connecting load transducer is secured to a 
Beome mOunting where an access hole through the beam housing 
aliows a set pin to lock or release the coupling for the body 
mounting. This decouples the moment indicating system when 
the test body's angle of rotation needs to be changed. This 
meen, reinstalled prior to testing, locks the test body to 
the front mount and thus allows twist to be transmitted through 
from the body. 

The rear end of the connecting gage is placed on the torque- 
beam mount and passes through a roller bearing support. This 
Support eeeeries the connecting load transducer to twist under 
the moment which may be experienced during the flow tests. 

The torque-beam mount (Figure 9) serves as the support and 
reterence position for the torque load gage, another 10 kg 
meree transducer which responds to deflections caused by the 
rotation of the rear end of the connecting gage. While the 
meomt Support of the torque transducer rests on the beam 
fount, the rear support is attached to an adjustable support 
Dracket. By shifting the connection position along the sup- 
port bracket, different orders of magnitude of moment can be 
compensated for and, therefore, permit more accurate recording 


mea calibration. 
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Prior to testing a new body, calibration of the drag and 
moment transducers was conducted. Known loads of up to 2,000 
meams were hung from the midlength of each body with the force 
Meeemin the drag position. The amplifier-recorder attenuation 
Me@tinegs were adjusted and compared for each different load to 
Smsoure consistency. 

Similarly, the moment transducer was calibrated by securing 
a clamp assembly, with a known moment arm, to the midpoint of 
the test body and then hanging various known loads from the arm 
to create test moment values. Calibration settings were checked 
Dy filling the tunnel with water and recording the signal gen- 
erated by the buoyant force. This signal was subsequently re- 
moved Dy zeroing the strip chart recorder before an actual flow 


test was conducted. 


SeeeoUSMARTNE-SHAPED TEST BODIES AND TESTING PROCEDURE 

Miree shapes, two cylinders and a flat plate, were tested 
during this research project. All of the bodies were 23.85 in. 
Seer Figures 10a and 10b provide cross-sectional profiles 
@eeecimensions for the 'D' and 'T' shapes, respectively. (Note 
that the letters 'D' and 'T' are used as designations only and 
do not necessarily define the shape of the cylinder cross- 
section.) Each of these cylinders was constructed from solid 
plexiglas, turned on a milling machine and then polished for 
mesmoOth finish. The length of each cylinder was cut so that 
a gap of approximately 0.06 in. was present between the tunnel 


wall and each end. 
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Peeeminum mounts were constructed on a lathe and then coun- 
Bemsunx and bolted into the cylinder ends so that the surfaces 
were smooth and flush. The attachment at the drag and lift 
end of each cylinder is round and compatible with the self- 
mepommne bearing. At the moment end, the cylinder attachment 
Momaetlat Cantilever beam which is compatible with the connec-. 
meemeco the front mounting of the connecting gage for the 
torque indicating assembly. 

The flat plate (Figure 10c) was fashioned from steel plate 
Memeo isned to a smooth surface. Its attachments are similar 
Mummmeee described above except for slots which hold the plate 
Meeeea to them. Each body is prevented from moving laterally 
Meepeeanecither wall of the tunnel by '0' rings attached between 
Mmemedantilever ends of the cylinder mounts and the force trans- 
Seer connections. (See Figures 6 and 8.) 

Beeeremoving the various locking devices previously de- 
memmped each test body can be rotated to various angles of 
Meetex, Both of the cylinders have a small pin hole at the 
mumeere an indicator fits for aligning the angle of rota- 
Meomeat 5 degree intervals. The indicator was lined up with 
the appropriate markings placed on the window surrounding the 
drag and lift gage. The indicator was removed from the cyl- 
Mee@er Shapes prior to testing. Flat plate angle of rotation 
memset DY using the knife-edge as a reference to the mark- 
ings on the window. 

Each body was tested at various angles of attack for drag, 


lift and moment measurement. At least three runs were conducted 
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Bemearag and lift, respectively, at each reference angle. 
Therefore, six moment runs were generated at each reference 
angle. Velocity was indicated by the elevation versus time 
muigtewhich was generated for every run. 

Raw analog data was produced from the force and moment 
meameaucers and the variable resistance probe by connecting 
their Outputs to individual Hewlett-Packard 7702B Recorder- 
Amplifiers and then to a Gould 2600 Strip-Chart Recorder which 
Was Operated at a speed of 50 mm per second for each run. 
Premres 11 and 12 are examples of the strip-chart recorder 
Weems) fOr the [-shaped cylinder and the flat plate, respec- 
tively. Note that for the purposes of presenting these ex- 
mmees, and the example in Figure 5, chart speed was reduced 
to 25 mm Pemescecoid, otrip €hart data was digitized using 
a Hewlett-Packard 9874A Digitizer coupled to an HP-9845A 
desktop computer. This analog data was processed in accord- 
ance with the governing equations described in Chapter IV, 
producing the graphical and tabular presentations of the 
BOrce-transfer coefficients included there and in the 


appendixes. 


D. FLOW VISUALIZATION 

Motion pictures, at film speeds of 32, 48 and 64 frames 
per second, were taken of the flow past the T-shaped body at 
mie zero reference position. The flow was filmed through the 
Window adjacent to the lift and drag gage housing. A plane 


of light across the cylinder was provided through slits 
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Figure 11. Analog Record for the T-Shaped Body at @ deg. 
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meeroximately one foot long cut into the front and rear walls 
Mempene tunnel. Shadow boxes were attached externally to the 
tunnel (Figure 1), into which were mounted high intensity 
imeants . 

After ensuring that the water in the tunnel was calm, 
polystyrene beads were dropped into the filled tunnel from 
Geemectly over the plane of the light through the cylinder 
cross section. After allowing the beads to sink down to the 
meeemeot the cylinder, while continuing to add more, the quick 
release valve was tripped to initiate the flow for the already 
running camera. 

me@eversiialization €Lforts were restricted to the T-shaped 
body because the view of the other shapes was obstructed by 


m@emadrag and lift gage housing and due to limited window area. 





ieee rosin RESULTS 


Pee PRESENTATION 

The experimental results (drag, lift and moment coeffi- 
cients [Cd, Cl, and Cm, respectively] for the D- and T-shaped 
Mmeemarical bodies; normal force coefficient [Cn] and the nor- 
malized moment arm [X/B] for the flat plate) are presented 
Setar ily in graphical form. Plots of the force-transfer 
coefficients versus an appropriate normalized displacement 
Mave been provided. For the purposes of discussion in this 
Mummers results are addressed, for a given body at a specific 
angle of attack, with reference to specific figures which dis- 
memercene plotted information. Since not all of the large 
vOlume of data are discussed herein, representative data 
together with the remainder of the graphs have been assembled 
in the appendixes by body shape and reference angle (angle of 
meeetek). (For more specific information refer to the List of 
PmeTidixes, p. 6). 

With a few exceptions, each plot shows data obtained in 


Pemee Separate experimental runs. 


Seer NITION OF THE FORCE-TRANSFER COEFFICIENTS 
he drag coefficient for the D- and T-shaped bodies is 


defined as 
ae » 
DRAG FORCE ( a) 
Cen 


Cd = 


on 





Mere € is the radius of the cylinder base circle; U, the 
velocity of the ambient flow; L, the length of the body; 
memo, the density of the fluid. Similarly, the lift coeffi- 
cient for these two bodies is defined as 


LIFT FORCE 
el Ue 


Cl (4) 


The moment coefficient 1s defined as 


MOMENT 


i) 
Ceo 


The positive directions of drag, lift and moment, along with 
M@empositive direction of angle of attack with respect to the 
Meeection Of the impulsive flow, are shown in Figures l5a and 
moomee lic normalized displacement for the D- and T-shaped 
bodies is UT/C. 

io Y Peeinit ion SeeEeoree and mement coefficients for the 
flat plate, a choice had to be made between seaprate lift and 
Meeeecociticients and a single, normal force coefficient. 
memececurrent literature on flow about flat plates often 
meeens tO the normal force and the use of one coefficient 
Meeeres Calculation of the lift or drag through the angle of 
attack, it was decided to present the normal force coefficient 
as 


FORCE NORMAL TO PLATE 
Bey 


oe (6) 


un 
IN) 





Lads 


DRAG (+) 


LIFT (+) 


Figure l3a. Definition Sketch for the D-Shaped Body. 
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Figure i3b. Definition Sketch for the T-Shaped Body. 





mere, 5 1S one-half of the plate width (see Figure l3c). For 
Gach angle of attack the lift and drag forces were separately 
Measured and vectorially combined to yield the normal force, 
Gueeae the Corresponding times. Only for 90 = 55 degrees, the 
memial tOrce Was measured directly, by suitably rotating the 
memce, transducer in order to compare the direct and indirect 
Seeetor Combination) methods of evaluating the normal force. 
Information regarding the moment acting on the plate is 

@eesented in terms of a normalized moment arm, X/B, which is 
@erined by 


X . MOMENT/Fn Gp 
B 


Ud| 


mumeriiy rn represents the distance, X, from the centerline 
mueemmor che plate to the point where the resultant normal 
force acts (see Figure 13c). The normalized displacement 
for the flat Dame ise 1/7 B. 

Note that "'angle of attack" is used synonomously with 
mimere Of rotation” as a reference position for the two cyl- 
andrical bodies. As shown in Figures 13a and 15b, this should 
not be confused with the angle of attack routinely associated 


Meme an airfoil. 


fee D-SHAPED CYLINDER 

mimrs bDOdy was tested at angles of attack of 0, +10, #20, 
meee and, -90 degrees. The drag, lift and moment coeffi- 
cients for the zero degree angle of attack are shown in Figures 


Be through 16, respectively, as a function of UT/C, the 
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Figure 1l3c. Definition Sketch for the Flat Plate. 
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Gememaeized displacement of the fluid. By carefully control- 
Mmemeecne flow mechanism in the vertical water tunnel, repeat- 
able and nearly constant velocities were maintained. The 


Beynolds number, Re = 2Uc/v, was, therefore, kept nearly con- 


it 


Stant at Re = 27830 cans, 

Mmematic effects are associated with the growth of the 
Meeeeenree or four vortices shed from the body following the 
Mermation Of the impulsive flow. At the start of the motion, 
Mm@emyortices grow rapidly and at a rate of growth such that 
Meemyorticity accumulates to an amount far in excess of that 
mumieein the later stages of motion. This excess vorticitv 
reduces the base pressure and causes a large drag overshoot. 
AS shown in Figure 14, the drag coefficient reaches a value 
Of slightly more than 1.6 with dual peaks in the region of 
Meeemoroot, extending from UT/C = 6 to 9. Following the shed- 
ding of the’second vortex (which can be estimated by studying 
muemeoot Of the lift coefficient, Figure 15), Cd drops rapidly 
and assympotically reaches a mean value close to that of a 
Me@emlar cylinder (in this case, Cd approaches a value of 0.8 
Mey C greater than 20). Although the experiments were not 
extended to UT/C values larger than about 26, the indications 
from the data in all cases are that beyond UT/C = 25 the drag 
coefficient will not change further beyond the assymptotic 
malues. 

For positive angles of attack the drag overshoot is also 


Meeoent with Cd reaching values of 1.7 to 1.8. At negative 
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amegles of attack, however, the drag coefficient increases 
meme r, to as high as 2.5 for 0 = -20 degrees. Note that 
@eemegative angles of attack, the flat face of the body is 
rotated toward the flow, presenting a sharp edge from which 
immediate separation can occur. The higher Cd values at 
negative angles of attack are, therefore, indicative of a 

more significant accumulation of excess vorticity in the early 
stages of motion (UT/C less than a value of about 11). (See 
memenaixes B through J). Assymptotic values of the drag coef- 
meetent for UT/C greater than 25 generally ranged from 0.8 to 
4. 

Peeure 15 shows the variation of the transverse force as 
Mepemct ion of UT/C. At the start of the fluid motion the lift 
mmee is directed toward the flat surface of the D-shaped body. 
(In this case Cl is negative by the definition adopted here 
and so is the lift force; see Figure 13a). Consequently, a 
Meemaped submarine with the flat surface horizontal (1.e., a 
level keel) is subjected to a negative lift force when pushed 
downward and a positive lift force when lifted toward the free 
Meeeeaee,. As displayed in Figure 15, the initial lift force 
resulting from the impulsive fluid motion is negative. A sub- 
Marine turning (to the right or left) would thus be subjected 
mena vertical force which tends to submerge the body for a 
lateral displacement of UT/C less than about 6. 

With reference to Figure 13a, the reasons for this behavior 


may be explained as follows. At the start of motion the flow 
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m_emeaces immediately at the sharp front edge of the flat 
Seerace. A vortex forms and grows rapidly. On the round 
meeenOt the D-shaped (bottom side of the submarine), however, 
M@emsecparation 1S not as immediate as on the top. Studies of 
Meemesively-started flows about circular cylinders have shown 
Meese paration does not start until the relative displacement 
Mumeeereaches a value of about 0.35 [4,62]. A circulation is 
imp<sed on the flow about the body which is of equal strength 
maemo posite sign for each shed vortex. Therefore, a vortex 
Meee ith a Clockwise circulation from the top of the D-shaped 
Meee Figure 13a) imposes a counter-clockwise circulation on 
Sie flow. As a result, the velocity increases and the pressure 
Meeredses On the round face relative to the other side where 
Muemerortex is first formed. The differential pressure creates 
a transverse force directed from the flat face (top) toward 
the center of the body. In other words, the lift force always 
acts on a body in a direction away from the last shed vortex 
and normal to the direction of the ambient flow. 

The subsequent shedding of vortices results in an alter- 
Meemae transverse force as shown in Figure 15. For the zero 
Meeree angle of attack, the lift force acting on the D-shaped 
Maye tiluctuates about a positive mean value. For all other 
angles of attack the value of the mean lift force is slightly 
Sueemecr tian that corresponding to the zero-degree angle of 


Beack (see Appendixes B through J). 
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ioe phenomenon of vortex shedding 1s non-stationary in the 
early stages of impulsive motion and consequently a single vor- 
tex shedding frequency does not exist. An estimate of the fre- 
quency was made from plots similar to Figure 15. Time intervals 
between consecutive peaks of the lift coefficient curves were 
Mereaeto derive a Strouhal number in terms of the period defined 
Dyetnose time intervals. Analysis of Figure 15 shows that the 
MmmeeetwO VOrtices shed from the D-shaped body, at zero degrees 
Meee Or attack, correspond to a Strouhal number, S = 2fC/U, of 
Mueeeecrmately 0.25. For other angles of attack, the Strouhal 
memeers ranged from 0.19 to 0.24 with no particular correlation 
Memeeite Specific direction or magnitude of the angle of attack. 

mmeeould not be inferred from the foregoing discussion 
Mmeeeine initial right or left turn of a D-shaped cylindrical 
body will necessarily and always result in a force which will 
MemaetO generate a diving motion. Other body- and control- 
Surface generated vortices will interact with the first vortex 
mooulting from tne impulsive motion of the submarine's turn. 
Meenough it was not the purpose of this investigation, anal- 
Memeomot the lift coefficient data presented herein indicates 
Beet the turn of a submarine may result in a diving or rising 
Memeo due in part to the impulsive cross-flow, the inter- 
action of all vortices in the flow surrounding the body, and 
Partly due to the heel of the submarine during the turn. 
Emphasis in this discussion has been placed on the separate 


effects of impulsive motion on the underwater Dody. 
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Figure 16 shows that the D-shaped body at zero degrees 
angle of attack is subjected to a positive moment as defined 
in Figure l3a. The direction of the moment could have been 
fererred from the vectorial sum of the lift and drag forces. 
Moment resulting from the impulsive flow tends to bank the 
Submarine in a direction opposite to the direction of heel 
meoeenormally occurs during the initial phase of the turn (a 
submarine rolls inboard during a turn [64]). In other words, 
Memeportion of the moment acting on the body which is created 
emeere impulsive cross-flow of the turn helps to stabilize the 
heeling position of the body. It is important to emphasize 
that the qualitative as well as quantitative understanding of 
the forces resulting from the impulsive motion is valuable in 
deriving a better overall understanding of all the forces 


maten determine the resultant motion of the body. 


Mel -oHAPED CYLINDER 

Meetees 17, 18 and 19 are the plots of drag, lift and 
moment coefficients, respectively, for the T-shaped body at 
Mero degrees angle of attack. Force and moment directions 
and other geometrical details are provided in Figure 15b. 
Meemkeynolds number, Re = 2UC/v, for the T-shaped cylinder 
was kept within a range of Re = 24380 + 240. 

AS in the case of the D-shaped body, the first vortex 
develops on the side where geometrical symmetry is distrubed, 
Meee, at the side of the protrusion. In the case of the D- 


shaped cylinder, the disturbance is the flat face (at the top 
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of the cylinder). On the T-shaped body geometric symmetry is 
meeeeanced by the presence of a flat-topped protrusion from 
mgemctop Of the cylinder (see Figure 13b). Rapid accumulation 
Meevorticity (referring to the case of the T-shaped body at 
mero degrees angle of attack) results ina large drag coeffi- 
Gient Over a relatively long time period or relative displace- 
meme cd = 2.2 at UT/C = 10.5). The terminal value of Cd in 
Figure 17 appears to approach a value of 1.5. Comparison of 
Saeures 1/7 and 18 shows that the shedding of the second vortex 
Seeiigsan abrupt reduction in the drag coefficient at UT/C = ll. 
The development of the third vortex from the protruding side 
Samune cylinder again increases the drag coefficient to a rela- 
Savely Smaller second maximum. Subsequent shedding of vortices 
may be expected to result in smaller fluctuations in the drag 
Meet icient. 

For increasing values of positive angles of attack (as 
@etined by Figure 13b), the assymptotic values of Cd are equal 
to about 1.9 until the angle is increased to +45 degrees where 
the assymptotic value (for UT/C greater than 25) drops to 1.4. 
meeceappendixes K through Y for the force-transfer coefficient 
Mmweaeat Other angles of attack). It is interesting to note 
mmmesene secondary maximum of the drag curve persists for all 
eemcaive angles of attack and that for positive angles greater 
than zero the drag force resulting from the development of the 
memrad vortex (at UT/C of about 14 to 15) is greater than the 


maximum reached during development of the first vortex. 
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Assymptotic values of the drag coefficient for negative 
mmeees OL attack up to -45 degrees are all in the range of 
Sempemerec. For © = -15 degrees and below, the large drag over- 
Mmmemeot the secondary maximum (third vortex development) per- 
ero as) it did for positive angles of rotation. This secondary 
fmmomum disappears for negative angles of attack greater than 
mmmewial to -20 degrees and, the drag coefficient character- 
meemes ClOSely approximate those of a circular cylinder. 

Peorezero degrees angle of attack, the initial lift force 
Seeeeetion is from the protrusion toward the axis of the cylinder 
Meet ive lift) during the development and shedding of the first 
Memetex.  inis 1S a characteristic which is similarly exhibited 
by the D-shaped cylinder. This phenomenon is clarified through 
m@emstuay Of the photographs depicting the evolution of the 
first two vortices (see Figures 20-26). Note that each photo- 
graph sequentially Gicmbays thestOrmatiton and shedding of the 
Mereemces fOr Corresponding relative displacements. <As seen in 
Figures 20 through 24, the first vortex, emanating from the 
protrusion on the T-shaped body, grows rapidly and then the 
Second vortex begins to grow on the opposite side of the cyl- 
maaer. Gradually the second vortex moves toward the down- 
memcam Stagnation point. The first 1ift maximum occurs at 
meemenil/C = S$, which corresponds to Figure 22. This is the 
Mistant at which the first vortex is still attached to (and 
about to be separated from) its shear layer and, the second 


memtex has not yet grown sufficiently. 
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Figure 28. 


Flow Rbout the T-Shaped Body at @ deg. 
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Figure 21. 


Flow About the T-Shaped Body at @ deg. 
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Figure 22. 


Flow About the T-Shaped Body 
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Figure 23. 


Flow Rboout the T-Shaped Body 
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Figure 24. 


Flow About the T-Shaped Body at dg deg. 
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Figure 25. Flow About the T-Shaped Body at 8 deg. and UTYC = 7.1. 
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Figure 26. Flow About the T-Shaped Body at 9 deg. and UTVC = 7.8, 





Meerigure 25, the first vortex begins to extend from the 
Mmear layer as the separation bubble of the second vortex 
(shed from the bottom of the body) has started to form. As 
time or the relative displacement of the fluid progresses, 
the second vortex grows and the lift coefficient passes through 
Mmmeeo at UT/C = 8, as seen in Figure 26 (where the first vor- 
tex is almost completely shed). At this point the circulation 
contributed by the developed vortices is nearly in equilibrium. 
Study of the subsequent frames of the motion pictures (not 
provided herein) has shown that the second vortex reaches full 
strength at UT/C = 10.3 and then the third vortex begins to 
develop from the protruding side of the body (where the first 
vortex was initiated). 

As discussed in detail in connection with the D-shaped 
Mepmmder at zero degrees angle of attack, the growth and motion 
of the first vortex sets DDE GiGcularton apout the body which 
Seeegual in magnitude and opposite in direction to that of the 
Mest vortex. This, in turn, results in a negative lift force 
(directed from the protrusion or the flat surface toward the 
axis of the body). The same result is evident for the T-shaped 
body at zero degrees angle of attack (see Figure 18) for the 
early stages of motion up to about UT/C = 10. 

Mes rOregoOing statements regarding the initial lift char- 
acteristics were true for all angles of rotation of the D- 
Shaped body. This does not imply that the same holds true 


for the T-shaped body. A detailed analysis of the lift 
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Me@eercient data presented herein provides a more complex lift 
meee action with exceptions to the general characteristics pre- 
viously discussed, which are dependent upon the angle of attack. 
Smerkxe the D-shaped cylinder, analog data from the force 
meemeaucer for the T-shaped body evidenced a brief but signif- 
Meme oeriod of dynamic change of lift direction within the 
Mamearate moments of the beginning of fluid motion. This 
Mumecde is indicated by the first few data points for UT/C less 
than 1 in Figure 18. Referring to Figure 13b and the remain- 
Mipeeiatt coefficient plots in the appendixes, this brief per- 
mere ron in the lift data involves that time when the flow is 
meemened and nearly irrotational. For positive angles of rota- 
meemeress than +10 degrees and through all negative angles of 
attack, the initial lift force at the very earliest instants 
Meeeuid motion is positive. In these cases, the attached 
meoweon the side of the protrusion must have larger velocities, 
nence lower pressures relative to the other side, and a small 
Meertave lift force of brief duration results. At positive 
angles of attack greater than +10 degrees, the initiation of 
fluid motion results in a small negative lift force for about 
Meenwsame relative displacement just described. Significant 
Meeeeton Of the protrusion past the vertical and closer to 
the approximate position of the downstream stagnation point 
apparently changes the initial velocity characteristics so 
that the higher velocity shifts to the bottom of the body and 
Phe brief duration of potential flow is characterized by nega- 


mive lift. 
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femee the presence of the protrusion on the body provides 
a consistent reference for the asymmetrical disturbance of the 
ambient flow, the direction and magnitude of the angle of 
meack Can also change the lift characteristics occuring during 
the period when the first vortex is shed (relative displace- 
Memes Of UT/C between 1 and 12). For all positive angles of 
Mmeede<, the lift force, relative to the shedding of the first 
vortex, begins with negative values (force directed from the 
Meeenwor the protrusion toward the cylinder axis), indicating 
Mumemoscparation on the side of the protrusion dominates the 
meeyestages of motion and the accumulation of vorticity. At 
Meeacive angles of rotation, however, where the protrusion is 
Meeacea tOward the flow, the reasons for strong, rapid genera- 
Moa the first vortex on that side appear to be lessened 
Mmemcne lift force changes direction. This is most evident 
as the negative angle of rotation is increased beyond 90 = -2Z0 
Meee (See Appendixes W, X and Y). A careful examination 
of the information provided by the plots of drag, lift and 
memernc coefficients for these angles gives some insight into 
the vossible reasons for this change. A relatively steep in- 
Meewoc Of the drag coefficient is evident, in each case, which 
1s directly related to the onset of asymmetry and separation. 
Mmewlitt force is negative for UT/C smaller than about 10 
during which time the first vortex grows at the side of the 
Meoerusion. When the T-shaped cylinder was rotated from -20 


©O -25 degrees angle of attack, the lift force became positive 
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Meeemescparation. The change in lift direction during the 
early stages of motion was more evident when the angle was 
Sememeased from -30 to -45 degrees. Without further investiga- 
tion through the use of flow visualization it is not possible 
memeontirm where separation first occurs. However, the effects 
of circulation applied to the flow about the body as a result 
Smee shed vortices, as previously discussed, imply that a 
clockwise circulation is present when the angle of attack is 
about -25 degrees. This would lead to the tentative conclu- 
Sion that the vortex shed from the bottom (round face) of the 
meayeis dOminant at these angles of rotation for values of 
Meetup to about 10 or 12. 

Two additional comments regarding the foregoing descrip- 
megumer the lift force characteristics of the T-shaped body 
with hespect to the angle of rotation are appropriate. 
mersc, the initial posmenem litt foree generated at the onset 
emeluid motion is very small and of such short duration that 
it would have negligible effects upon the motion of an under- 
Water body, with this cross-section, undergoing an impulsive 
Siem. the discussion regarding this particular experimental 
observation is presented here because it provides further ex- 
perimental information relevant to the applications of poten- 
meal flow theory. Second, the explanation of the shift in 
Mmerection of the initial lift force as a function of increas- 
ing negative angles of attack is somewhat conjectural and re- 
quires further experimental investigation, with detailed flow 
Sesualization, for its substantiation. 
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Mabike the D-shaped body, the shedding of the vortices 
Causes significant changes in the moment acting on the T- 
Smempea DOdy (see Figure 19). Tnis is due primarily to the 
fact that the protrusion increases the available moment arm 
whereas, with the D-shaped body, the flat face decreases the 
Moment arm. Consequently, underwater bodies with protrusions 
meeemore likely to undergo larger amplitudes of moment fluc- 
tuations as a result of an impulsive motion. This may give 
Pise tO greater stability problems. Representative tabulated 
mmeamand plots of the force-transfer coefficients are pre- 
sented in Appendixes K through Y. 

Vortex shedding frequencies and the corresponding Strouhal 
numbers for the T-shaped body have been examined in a manner 
Similar to that for the D-shaped body. [It was found that S = 
Meeereror zero degrees angle of attack. The Strouhal number 
for all other CAG komo rOratren varied from 0.15 to 0.21. 
Decreasing Strounal numbers, relative to the D-shaped body, 
are in conformity with the general rule that the Strouhal 
number decreases with the increasing degree of bluffness of 
the body. 

A summary observation is in order regarding the descrip- 
mom Of the characteristics of the forces acting on and the 
vortices shedding from the D-shaped and T-shaped bodies. 
First, the asymmetry of the body relative to the direction 
of the flow removes the randomness from the position of the 


Mreowtn Of the first vortex. In the case of a circular 
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Meeemaer the first vortex may form from either side, at random. 
The introduction of a geometric deformity to the otherwise 
memenmmetrical shape of the cylinder (such as a flat face or 
meemotrusion) fixes the direction of the lift force and the 
meeeeton of the first vortex, in the early stages of the 
motion (as evidenced by the data and their repeatability for 
multiple runs). Second, an enlargement rather than a reduc- 
tion of the otherwise circular cross section increases the 
moment and accentuates the fluctuations of the moment due to 
the vortex shedding. Furthermore, an increase in the bluff- 
Mess Of the body, by the addition of a protrusion, decreases 
the Strouhal number and hence the vortex shedding frequency. 
Mmeeery, tne interaction of the body and the vortices in the 
Smrly Stages of the impulsive motion is such that a satisfac- 
memye analytical model requires both the representation of the 
shear layers by discrete vortices and, the evolution of the 
meme-dependent boundary layer. Future analytical efforts 
Seema find the data presented herein useful for their evolu- 
meomeand Calibration. Although the results presented here 
Were obtained at nearly constant and relatively low sub- 
Critical Reynolds numbers, the observations and conclusions 
Meeenot expected to differ significantly from those cited 
above even at much higher subcritical Reynolds numbers. One 
must add that the upper limit of the supercritical range is 
meee Known fOr non-circular bodies and can be determined only 


merough painstakingly difficult experiments. 
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fee LAT PLATE 

ieee geometrical characteristics of the flat plate and the 
Setection of the forces and moment acting on it are shown in 
Meeere 15c. The Reynolds number for this body was kept at a 
Meaniy constant value of Re = 2UB/v = 20750 + 260. The results 
m@eenwaiscussed with respect to two representative angles of 
attack (9 = 90 and 55 degrees), in order to accentuate the 
differences brought about by the shedding of the vortices on 
M@eemnonrmal force coefficient (Cn) and the position of the re- 
Sultant force as defined by the normalized moment arm (X/B). 

Figure 2/7 shows the normal force coefficient data for 90 = 
90 degrees. As anticipated, the symmetric growth of the vor- 
tices causes a large normal force coefficient to develop dur- 
Mmeethe early stages of the fluid motion. For UT/B larger 
than about 12, Cn reaches a nearly constant value of about 
oe. With the plate positioned normal to the ambient flow, 
mgessnedding of the vortices does not cause large fluctua- 
mons in Cn. 

Figures 28 and 29 are the plots of the normal force coef- 
ficient for 90 = 55 degrees. The results shown in Figure 28 
Mere obtained by summing vectorially the lift and drag forces. 
Data plotted in Figure 29 were obtained by measuring the 
normal force directly, as described in Chapter 3. A quick 
review of the two figures shows that they are nearly identical. 
Thus the indirect method (vector sum) of evaluating Cn is free 


from experimental errors within acceptable limits. This fact 
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is further evidenced by Figure 30 where the results of two 
mepresentative runs for each method are superimposed. 

maeemring £O Figure 30, at the early stages of fluid motion 
Cn increases dramatically to values never before noted in the 
relevant literature. The shedding of the first vortex at 
memmtl/B = 9 (also refer to Figure 31) brings about an 
memnoe Change in Cn. A subsequent maximum in the normal 
merece coefficient occurs at UT/B = 14, as a result of the 
Mmmem@eine of the second vortex. In fact, a flat plate at all 
angles of attack, other than those where the plate is nearly 
normal to the flow (90 = 80 and 85 degrees, see Appendixes 2 
and AA) exhibited the largest Cn value at about UT/B = 14. 

Wathin the range of relative displacement encountered 
during these tests, Cn exhibits smaller fluctuations and 
reaches nearly constant values for relative displacements 
greater than UT/B = 14. The results of discrete vortex 
analysis [54] for 9 = 50 degrees, show normal force coeffi- 
Memes With a maXimum of 2.5. Constrasted with a value of 
Meeencountered at the corresponding angle of attack for 
these experiments (see Appendix AD}, the present experimental 
mamaence implies that the predictions of the discrete vortex 
model are in need of further fine-tuning and correction. 

The normalized moment arm (X/B) for an angle of attack of 
memcdegrees 1S Shown in Figure 31. Cleariy the fluctuations of 


he moment arm are directly related to the shedding of the 


Mertices. Generally, the maximum absolute change in the 
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memenc arm decreases with a decrease in the angle of attack. 
Mmemecne plate 1s nearly normal to the flow (0 = 85 degrees 
Or more) the moment arm is reduced and remains closer to the 
Mumeeecenterline. As evidenced by the data presented in 
Mmeemaix 2, at very large angles of attack, the normalized 
moment arm and the normal force coefficient are noticeably 
Mememscensitive to changes in the ambient flow conditions 
Meemeat Smaller angles of attack. When the angle of attack 
is reduced to less than 40 degrees the normalized moment 
Mameremains positive throughout the range of relative dis- 
megeement encountered herein and there is a significant 
Meomction in the fluctuations of X/B in response Pomme 
memeex-shedding. 

The analog traces of the lift force and moment and the 
plots of the normalized moment arm were sufficient to identify 
the vortex shedding times and to estimate the Strouhal number, 
Mueeete/U, for most angles of attack. For 9 = 55 degrees 
@emeure 31) the Strouhal number was about 0.22. The Strounhal 
mumper ranged from 0.14 at 9 = 85 degrees to 0.25 at © = 40 
Meerecs. Clearly, the Strouhal number decreases as the bluff- 
ness of the body is increased, by increasing the angle of 
merack. it is also important to note that the fluctuations 
in moment acting on the flat plate are much more dramatic 
Mean for cither the D- or T-shaped bodies. This is partly 
Mesto the fact that the vortex strength for a flat plate 


Memerareecr and partly oécause the sharp edges of tne plate 
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fix the separation points and increase the spanwise coherence 
Mmeemerene plate. For bodies without sharp edges, the separa- 
meen points are mobile and there are often significant phase 
Mees in pressure at the corresvonding points along the 
plate. This leads to a reduction in the spanwise coherence 


Mmemeene vortices and hence to a reduction in lift forces. 


Meer eCAPITULATORY SYNTHESIS 

Bae hydrodynamical behavior of the three bodies under con- 
Sideration may be compared relative to each other, relative 
femeteat Of a circular cylinder and, relative to analytical 
eemumerical results. 

An extensive literature search and the perusal of the 
most relevant works have shown that the analytical or numer- 
Heal methods are not yet in a position to describe the behav- 
mepeor the large scale motion in the wake of bluff bodies 
Subjected impulsively started flow. Solutions based on the 
Navier-Stokes equations are limited to extremely small 
Peynolds numbers and to relatively small displacements. 
Solutions based on the discrete vortex method, or on its 
Variations, suffer from the problems associated with the 
Meeroduction of the nascent vortices, vortex excursions, 
Shear-layer instabilities, inaccurate representation of 
Mescous and turbulent dissipation of vorticity, and numer- 
Ous other numerical instabilities. The assumption of an 


inviscid fluid imposes an unrealistically large number of 
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Bemareions on the simulation of the behavior of a viscous flow. 
{t does not, therefore, appear that the numerical methods 
mepeeted Cither directly or indirectly (through the use of one 
or more conformal transformations) are in a position to pre- 
Mmeemcne fOrce, moment and vortex-shedding characteristics of 
meatier bodies in time-dependent flows even for the most manage- 
able impulsively-started flow. Experiments must be carried 
Out to delineate the most important features of such flows 
partly for practical purposes and partly, and perhaps more 
memorcantly, for the purpose of providing physical insight 
toward the development of analytical metnods. 

Of the three bodies, the T-shaped body had mobile separa- 
tion points, the flat plate had clearly fixed separation 
points, and the D-shaped body had at times only mobile separa- 
meen points and at times both mobile and fixed separations, 
@epending on the angle of attack. In any case, however, the 
D-shaped cylinder had at least two and as much as three sep- 
Merenon points. It was for the purposes of accentuating the 
Beeeects Of these three types of separations that these three 
particular body shapes were chosen. 

The results have shown that the fixing of the separation 
point has numerous and yet incalculable consequences. The 
vortex strengths are larger (as evidenced by the magnitude 
meecm@e rorce fluctuations), the spanwise coherence is in- 
creased out of necessity, and the Strouhal number is reduced 


partly because of the increase of the strength of the vortices 
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and partly because of the bluffness of the body. On the other 
hand, a body with mobile separation points, with a preferred 
meertion Of initial separation (T-shaped cylinder), exhibits 
Mmeesacavely smaller drag overshoots and fluctuations of forces 
and moment. The magnitude of the moment obviously depends on 
Miewamount of protrusion relative to the base circle of the 
memiaer Shapes. the protrusion also fixes the direction of 
the lift force in the initial stages of the motion. This then 
as the basic difference between the impulsive flow on a circu- 
meeecylinder and on a T-shaped cylinder, both of which have 
[Meeule Separation points. On a circular cylinder the side 
from which the first vortex sheds is random and determined 
@emt@arily by small disturbances in the ambient flow. Fur- 
Mumemere, the circular cylinder experiences little or no 
moment. On a geometrically deformed cross section, as in 
i@emcase Of the T-shaped body, the randomness is removed 
meeom the shedding of the first vortex, the direction of the 
Merce is tixed and a moment is created. At large angles of 
Meee (€.¢g., -25 degrees) the T-shaped body may exhibit 
three mobile separation points. This may lead to a change 
imieche direction of the initial lift force, but still does 
Meee introduce randomness into the evolution of the lift 
morce. 

her a Dbody with mixed separation points, as in the case 
Ot the D-shaped cylinder, the angle of attack plays a rela- 


tively more important role in determining the flow characteristics. 
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Sec) again the position of the first vortex and the direction 
Seemte initial lift force are fixed on the side of the body 
where the separation 1s fixed. A closer examination of this 
and other bodies shows that any geometrical deformity which 
mmeeripitates the occurrence Of separation necessarily leads 
memerepreferential shedding of the first vortex and to the 
Seememmination Of all other features of the flow. On a cir- 
Gular cylinder the separation does not start immediately, 

but at a fixed pertubation in the otherwise smooth circular 
cross-section, separation of the flow takes place immediately. 
Consequently, any physical disturbance imposed on the body, 
De it a sharp edge, a rounded protrusion or a flat surface, 
Meeetpitates separation either relative to a circular cylinder 
Smeerelative to the undisturbed side of the body, and fixes 
the rest of the history of the motion by removing any degree 


of randomness. 
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VeeeecCONCEUS IONS 


The experimental investigations of the impulsively-started 
flow about a D-shaped cylinder, a T-shaped cylinder anda flat 
eee resulted in the following conclusions: 

1. Shedding of the first two or three vortices from the 
Beeama I-shaped cylinders produced a significant drag over- 
Shoot of as much as fifty percent. 

2. Drag coefficients for the D-shaped cylinder assymptoti- 
Gally approached values in the range of 0.8 to 1.4 for UT/C 
Bmeacer than about 20. For the T-shaped body, assymptotic 
Values of Cd ranged from 1.4 to 1.9. There were no signif- 
meate tluctuations in the drag coefficients for these bodies 
Merelative displacements beyond UT/C = 20. 

Eemeiitcher Cd values occur, for the D-shaped body, at 
negative angles of rotation. 

Geerconuetne D-shaped cylinder the lift force direction is 
meeeterve during the early stages of motion, for all angles of 
attack investigated here. 

ceueeor the T-shaped body, at most angles of attack, Cd 
reached a secondary maximum due to the development and subse- 
Guent shedding of the third vortex. This secondary maximum 
Mumerreer than that occtirring during the earlier stages of 


the motion. 
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Meee lne Lift force direction during the early stages of 
impulsive motion about the T-shaped body is dependent upon 
the angle of rotation relative to the ambient flow. 

7. Asymmetry of a body relative to the direction of the 
impulsive flow removes the randomness from the position of 
@memeration of the first vortex and fixes the initial direction 
Meede Lift force acting on the body. 

feeeeosymmetric enlargement of an otherwise circular cross- 
meerron increases the magnitude of the moment experienced by 
tne body and accentuates the fluctuations in that moment. 
Moment fluctuations for the flat plate are more dramatic 
Bem £Or eC€ither of the cylindrical bodies. 

9. The maximum normal force coefficients for the flat 
Meete are Significantly greater in magnitude, and in the ampli- 
mmeenot fluctuations, than those predicted by the discrete 
vortex model. 

Mmemeocrounal numbers ranged from 0.19 to 0.24 for the D- 
meme cylinder, from 0.13 to 0.21 for the T-shaped cylinder, 
mmemetrom 0.14 to 0.25 for the flat plate. 

11. There are, at present, no reliable analytical techniques 
which accurately predict the behavior of the impulsive flow 
Mout bluff bodies. Those available have proven to be defi- 
meme in one Or more ways. Experimental measurements and 
observations still form the basis for a comprehensive under- 
Standing of the behavior of impulsively-started flow about 


mitt bodies. 
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VI. RECOMMENDATIONS FOR FURTHER RESEARCH 


It is recommended that additional research on impulsively- 
Started flow about bluff bodies be conducted in the following 
aac a : 

Mupeempulsive flow tests of the three bodies presented 
herein should be conducted at significantly different, and 
meecerabDly higher, Reynolds numbers. 

meeactempts should be made to determine the velocities 
and vortex strength either through visualization or through 
meemuse Of a laser velocimeter. 

Meee the pressure distribution about the circumference of 
the bodies should be measured and evaluated as a function of 
memes trom the start of the fluid motion. 

@m—mexecensive flow visualization photographs of the tests 
on these bodies should be conducted with attempts to monitor 
Meeectow for UT/C values of 100 or more. Accurate correlation 
Meme time-sequence of the photographs with the analog traces 
Meetne torces and moments should be an important part of this 
Setort. 

5. The present results should be extended to larger UT/C 
values (in connection with the recommendation above) in order 
Meeverity the assymptotic behavior of the drag coefficients 


Beato Calculate the Strouhal numbers for larger times. 


O7 





Seine present study should be extended to different 
Mross-sections. Rectangular cylinders present the possi- 
feet y Of multiple separations and reattachments of the flow. 

7. The T-shaped body should be tested with its axis at 
different angles of pitch to determine the effects of the 
body vortices and the interaction of all vortices shed when 
the body's axis is inclined other than normal to the flow. 

8. hese investigations should be expanded to cover gen- 
eral three-dimensional bluff bodies, on which there is very 
Meeele Or no research, save for the sphere. 

Oo. The effect of ambient flow turbulence on the impul- 
Mme -Started flow should be investigated in order to under- 
Stand the initiation of vortex asymmetry on axisymmetric bodies. 

10. The foregoing investigations should be extended to 
mematively more complex time-dependent flows where the velocity 
varies as a known function of time or displacement in order to 
mm@@erstand the effect of the local acceleration of the ambient 
Sroweon the vortex development and on the time-dependent 


resistance. 
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RPPENDIX AR: REPRESENTATIVE DATA FOR FLAT PLATE AT 88 AND 75 DEGREES 
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APPENDIX AD: REPRESENTATIVE DATA FOR FLAT PLATE AT 58 AND 45 DEGREES 
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RPPENDIX RE: REPRESENTATIVE DATA FOR FLAT PLATE AT 46 AND 38 DEGREES 
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